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Researchers at the National Institute of Standards and Technology have been using a watt balance, 

NIST-3, to measure the Planck constant h for over ten years. Two recently published values disagree 
by more than one standard uncertainty. The motivation for the present manuscript is twofold. 

First, we correct the latest published number to take into account a recently discovered systematic 
error in mass dissemination at the Bureau International des Poids et Mesures (BIPM). Second, we 
provide guidance on how to combine the two numbers into one final result. In order to adequately 
reflect the discrepancy, we added an additional systematic uncertainty to the published uncertainty 
budgets. The final value of h measured with NIST-3 is h = 6.626 069 36(37) x 10 -34 J s. This 
result is 77(57) x 10 -9 fractionally higher than hgo. Each number in parentheses gives the value of 
the standard uncertainty in the last two digits of the respective value and hgo is the conventional 
value of the Planck constant given by hgo = 4/(7\'j_ 90 i?K-9o), where Tfj-go and Rk -90 denote the 
conventional values of the Josephson and von Klitzing constants, respectively. 


I. INTRODUCTION 

Researchers at the National Institute of Standards and 
Technology (NIST) have published two results for the 
Planck constant, ho? in 2007— and /ii 4 in 20142 using the 
same apparatus, the third generation of the NIST watt 
balance, referred to as NIST-3. The published values are 

h 07 = 6.626 068 91(24) x 10“ 34 J s, or 
^Z. - 1 = 8(36) x 10" 9 , 

h 90 

and 

h 14: = 6.626 069 79(30) x 10" 34 J s, or 

-1 = 141(45) x 10" 9 , 
hg o 

respectively. Each number in parentheses gives the value 
of the standard uncertainty (k = 1) in the last two 
digits of the quoted value. The symbol h 90 denotes 
the conventional value of the Planck constant given by 
hgo = 4/(A" 2 _ 90 R K -9o), where A'j -90 and Ak -90 denote 
the conventional values of the Josephson and von Klitz¬ 
ing constants, respectively. 

Before we discuss the two measurements, we note that 
the 2014 result must be adjusted due to an offset in the 
SI unit of mass disseminated by the Bureau International 
des Poids et Mesures (BIPM). During the extraordinary 
comparison^ of the international prototype of the kilo¬ 
gram (IPK) at the BIPM, it was found that the mass 
unit as maintained by the BIPM is 35 pg larger than the 
mass of IPK—. 

All measurements of the Planck constant discussed 
here were performed with the NIST platinum-iridium 
prototype No. 85, known as K85. Figure [T| shows the 
calibration history of K85. The first calibration certifi¬ 
cate for K85 was issued by the BIPM in November 2003. 


Shortly after being put in service with NIST-3, the mass 
of this prototype drifted upwards with a rate of approx¬ 
imately 5/ig per year. At the time of the NIST-3 exper¬ 
iment, this drift was taken into account in the calcula¬ 
tions of the published values of h. In 2010, the United 
States shifted its mass scale by 45 pg/kg, because it was 
found during a routine calibration of the US national 
standard prototype K20 that the BIPM mass scale and 
the US mass scale differed by 45 pg/kg— At the end of 
2011, K85 was sent for calibration to the BIPM, where 
it was washed two times. The measurements of the mass 
of K85 performed by the NIST mass and force group 
and the BIPM in 2012 and 2013 are in agreement within 
uncertainties. A constant value for the mass of K85 of 
1 kg — 738.3 pg was used to calibrate all NIST-3 measure¬ 
ments made in 2012 and 2013. 

In 2014, during the extraordinary mass comparison, it 
was discovered that the mass unit as maintained by the 
BIPM is 35 pg/kg larger than the mass scale set by the 
IPK—. It is believed that the difference between the mass 
unit as maintained by the BIPM and the IPK has accu¬ 
mulated between 2000 and 2014. Calibrations made at 
the BIPM before 2000 are therefore deemed good cali¬ 
brations with respect to the IPK—. But, calibration cer¬ 
tificates issued by the BIPM after 2000 have to be cor¬ 
rected to be consistent with the international definition of 
the kilogram. According to a numerical model made by 
the BIPM, the calibration value assigned to K85 by the 
BIPM in 2003 must be decreased by approximately 4 pg 
and the one assigned in 2012 by approximately 35 pg— 
The latter correction largely erases the adjustment made 
to the US mass scale in 2010. 

After K85 had arrived at NIST at the end of 2003, it 
was measured against the national prototype, K20, by 
the NIST mass and force group. The national prototype, 
as well as the national check standard, K4, have been 
calibrated at the BIPM in 1999. At this time the mass 
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unit as maintained by the BIPM was still synchronized to 
the IPK. Hence, it is believed that the calibration values 
determined by the NIST mass and force group for K85 
from 2004 through the end of 2011 are accurate repre¬ 
sentations of the mass of K85 with respect to the mass 
of the IPK. The h measurements taken during this time 
do not need to be corrected^. 

In 2011, the mass and force group shifted its mass scale 
to be consistent with the mass scale given by the mass 
unit as maintained by the BIPM, which, as we now know, 
was different than the mass unit given by the IPK. Hence, 
all data taken with NIST-3 in 2012 and 2013 used a value 
for K85 that was tied to the mass scale as maintained by 
the BIPM, which was about 35 pg different from the 
mass of the IPK at that time. Therefore, the measured 
value for hu needs to be corrected down relatively by 
35 x 10" 9 . This yields 

/i 14corr = 6.626 069 56(30) x 10 “ 34 J s, or 
- 1 = 106(45) x 10” 9 . 

«90 

In summary, the time dependent shape of the correc¬ 
tion to the NIST h values is a step function rather than a 
smooth interpolation. The relative correction that must 
be applied to the published h values is 0 x 10~ 9 for data 
taken before Dec. 2011 and —35 x 10 -9 for data taken 
thereafter. 
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FIG. 1. The mass of the platinum-iridium prototype No. 85 
(K85). The circles represent calibrations that were performed 
at the time of the experiment, either by the BIPM mass de¬ 
partment (solid) or the NIST mass and force group (open). 
The diagonal crosses represent values that were calculated by 
applying a correction issued by the BIPM^ to the original 
BIPM measurements. The upright crosses denote values that 
were calculated from measurements obtained by the NIST 
mass and force group by subtracting a correction. The correc¬ 
tion is 0 pg before December 2010 and 35 pg after December 
2010. The upright crosses represent the current best estimate 
of the mass of K85 relative to the IPK. 


January 2012 are corrected by the 35 x 10 ~ 9 as discussed 
above. The lower panel shows the data grouped in daily 
data runs. The data can be divided into three epochs. 
The upper panel shows the values obtained by taking 
the mean of the data in each epoch. The error bars indi¬ 
cate the statistical standard deviation of the data. The 
mean value and the statistical standard deviation for each 
epoch are also given in table Q] 


-c 


-c 


150 
100 
50 
0 

200 
150 
100 
50 
0 
-50 

'03 '04 ’05 '06 '07 '08 '09 '10 '11 '12 '13 '14 
Time (year) 


Epoch 1. Epoch 2 Epoch 3 



FIG. 2. Bottom Panel: All measurements obtained with 
NIST-3 using the platinum-iridium prototype No. 85. Top 
Panel: The average value obtained during each epoch. The 
error bars give the statistical standard deviation of the data. 
The data points in the third epoch in both panels are cor¬ 
rected down relatively by 35 x 10~ 9 from the original data 
to account for the offset between the IPK and the mass unit 
as maintained by BIPM that was discovered at the extraordi¬ 
nary comparison. The solid black lines in both panels give the 
unweighted average, 77 x 10~ 9 , of the values from the three 
epochs listed in table 1. The gray bands around the solid lines 
indicate the one sigma uncertainty of the final value reported 
here. 


The average value of the data in the first epoch is 
slightly larger than the value published ini. This is due 
to the fact thati includes data up to July 2006, while 
the first epoch includes data up to the end of 2009. The 
additional data increased the combined mean by a non¬ 
significant amount, i.e., the mean of the epoch and the 
published value are consistent within uncertainties. 


epoch time ( 7^0 — l) x 10 9 7T x 10 9 


1 2004-2009 

29 

19 

2 2010/2011 

95 

37 

3 2012/2013 

106 

38 

unweighted mean 

77 


TABLE I. The average values of ( h/hgo — 1) x 10 9 for the three 
epochs. Here ah denotes the statistical standard deviation in 


each epoch. 


The relative difference between the values of /107 and 
/ii 4 corr is 98 x 10 -9 . All data obtained with NIST-3 us¬ 
ing K85 are shown in Figure [2] The data points after 


The measured value of h shifted by a relative amount 
of approximately 70 x 10~ 9 , which at the time was not 
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recognized, but we now term as the beginning of the sec¬ 
ond epoch of data recording. Nearly simultaneously with 
the (only later realized) long term shift, the day to day 
scatter of the data increased by a factor of two. At the 
time, the experiment was being incrementally improved 
while being studied to determine the cause of the sud¬ 
den changes. Even recent reexamination of the change 
records and effects could find no clear, single cause for 
the shift or the increase in scatter. 

Since 2010, the measured values were stable around a 
mean value of h/h 90 — 1 ~ 100 x 10 -9 and the data could 
have been summarized into one longer epoch. However, 
significant changes were implemented starting in 2012 as 
summarized inS.. Most importantly, the 2012 and 2013 
measurements were performed blind with the researchers 
involved not knowing the true result until the unknown 
offset was publicly revealed in June 2013. Considering 
the fact that major changes between the second and third 
epoch did not alter the mean values gives us more con¬ 
fidence in the later number. Moreover, the shift in the 
mass scale discussed above can be interpreted as a second 
blind bias added to the experiment. Yet, the 2012 and 
2013 results agree with the measurements taken between 
2010 and 2011. Hence, we decided to give the later data 
more weight (two epochs after 2010 vs. one epoch be¬ 
fore 2010). Weighting the data after 2010 more heavily, 
also reflects the fact that our understanding of the ap¬ 
paratus grew as a function of time. Towards the end of 
the data shown in Fig. [2j we were more confident in our 
understanding of the system. 

After more than fifteen years of experience with this 
apparatus, we believe that our best measurement of h is 
obtained from the unweighted average of the values ob¬ 
tained in the three epochs listed in Table 1. This average 
is h/h 90 - 1 = 77 x 10 - 9 . 

Despite all our experience, we acknowledge that we do 
not understand the cause of the approximately 70 x 1(W 9 


relative shift. This lack of understanding must be re¬ 
flected in the uncertainty assigned to the final value. We 
assign half of the observed shift, 35 x 10 -9 , as an ad¬ 
ditional relative uncertainty component to account for 
a possible unexplained systematic effect in the measure¬ 
ments. By adding this type B component in quadrature 
to the relative uncertainty published in^, a combined rel¬ 
ative standard uncertainty of 57 x 10~ 9 is obtained. 

The final value for the Planck constant obtained with 
NIST-3 in ten years of measurement is 

hmsT-s = 6.626 069 36(37) x 10~ 34 J s, or 
/tNIST ' 3 - 1 = 77(57) x 10 -9 . 

“90 

A previous version of the NIST watt balance (NIST-2) 
was used to determine a value of the Planck constant in 
a measurement campaign lasting four months in 1998^. 
The result 

/inist -2 = 6.626 068 91(58) x 10" 34 Js, or 
/lNIST ~ 1 2 3 4 5 - 1 = 8(87) x 1CT 9 

“90 

was obtained using two gold masses with a combined 
mass of 1 kg. Hence, the value /inist -2 is independent 
of the mass changes discussed above and there is no need 
to modify this value or its uncertainty. Based on the 
previous estimate^ of 0.14 for the correlation coefficient 
between /inist -2 and we estimate 

r(h N iST-3> /inist-2) = 0.09 

to be an upper limit for the correlation coefficient be¬ 
tween /inist -3 and /inist- 2 - 

We thank the mass and force group at NIST, especially 
Zeina Kubarych and Patrick Abbott, for calibrating our 
masses and for sharing their expertise in mass metrology 
with us. 
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